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Abstract 

The type la supernovae observational data are used to estimate the parameters of a cosmological 
model with cold dark matter and the Chaplygin gas. This exotic gas, which is characterized by 
a negative pressure varying with the inverse of density, represents in this model the dark energy 
responsible for the acceleration of the Universe. The Chaplygin gas model depends essentially on four 
parameters: the Hubble constant, the velocity of the sound of the Chaplygin gas, the curvature of 
the Universe and the fraction density of the Chaplygin gas and the cold dark matter. The Bayesian 
parameter estimation yields Ho = 62.1±l% km/Mpc.s, £l k0 = -0.84t{;ll, fi m0 = O.Oj;^® 2 , Q c0 = 
1.40^} 5 i6, A = c 2 = 0.93l°;!£c, to = u'.2±f%Gy and q = -0.98t.lH. These and other results 
indicate that a Universe completely dominated by the Chaplygin gas is favoured, what reinforces 
the idea that the Chaplygin gas may unify the description for dark matter and dark energy, at least 
as the type la supernovae data are concerned. A closed and accelerating Universe is also favoured. 
The Bayesian statistics indicates that the Chaplygin gas model is more likely than the standard 
cosmological constant (ACDM) model at 55.3% confidence level when an integration on all free 
parameters is performed. Assuming the spatially flat curvature, this percentage mounts to 65.3%. 
On the other hand, if the density of dark matter is fixed at zero value, the Chaplygin gas model 
becomes more preferred than the KCDM model at 91.8% confidence level. Finally, the hypothesis 
of flat Universe and baryonic matter (fi&o = 0.04) implies a Chaplygin gas model preferred over the 
KCDM at a confidence level of 99.4%. 



PACS number(s): 98.80.Bp, 98.65.Dx 



1 Introduction 

The combined data from the measurements of the spectrum of anisotropies of the cosmic microwave 
background radiation pQ and from the observations of high redshift type la supernovae Ej indicate 
that the matter content of the Universe today may be very probably described by cold dark matter 
and dark energy, in a proportion such that £ld m ~ 0.3 and Q^e ~ 0.7. These dark components of the 
matter content of the Universe manifest themselves only through their gravitational effects. At same 
time, a fraction of the dark components agglomerates at small scales (cold dark matter) while the other 
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fraction seems to be a smooth component (dark energy). The dark energy must exhibit negative pressure, 
since it would be the responsible for the present acceleration of the Universe as indicated by the type la 
supernovae observations, while the cold matter must have zero (or almost zero) pressure, in order that it 
can gravitationally collapse at small scales. 

The nature of these mysterious matter components of the Universe is object of many debates. The cold 
dark matter may be, for example, axions which result from the symmetry breaking process of Grand Uni- 
fied Theories in the very early Universe. But, since Grand Unified Theories, and their supersymmetrical 
versions, remain a theoretical proposal, the nature of cold dark matter is yet an open issue. 

A cosmological constant is, in principle, the most natural candidate to describe the dark energy. It 
contributes with an homogeneous, constant energy density, its fluctuation being strictly zero. However, 
if the origin of the cosmological constant is the vacuum energy, there is a discrepancy of about 120 orders 
of magnitude between its theoretical value and the observed value of dark energy ■ This situation can 
be ameliorate, but not solved, if supersymmetry is taken into account. Another candidate to represent 
dark energy is quintessence, which considers a self-interacting scalar field, which interpolates a radiative 
era and a vacuum dominated era |SJ |7| . But the quintessence program suffers from fine tuning of 
microphysical parameters [Sj. 

Recently, an alternative to both the cosmological constant and to quintessence to describe dark energy 
has been proposed: the Chaplygin gas O^H^DEl The Chaplygin gas is characterized by the equation 
of state 



where A is a constant. Hence, the pressure is negative while the sound velocity is positive, avoiding 
instability problems at small scales |13|. The Chaplygin gas has been firstly conceived in studies of 
adiabatic fluids |14| . but recently it has been identified an interesting connection with string theories 
[151 1161 IT7| . Some extensions of the Chaplygin gas have been proposed ^2] through the equation of state 



However, the connection with string theories is lost unless a = 1, and for this reason we will only consider 
in what follows the equation of state Q. 

Considering a relativistic fluid with the equation of state , the equation for the energy- momentum 
conservation relations leads, in the case of an homogeneous and isotropic Universe, to the following 
relation between the fluid density and the scale factor a: 



where B is an integration constant. This relation shows that initially the Chaplygin gas behaves as a 
pressureless fluid, acquiring later a behaviour similar to a cosmological constant. So, it interpolates a 
decelerated phase of expansion to an accelerated one, in way close to that of the quintessence program. 

In this work, we will constrain the parameters associated with the Chaplygin gas model (CGM) using 
the type la supernovae (SNe la) observational data. Specifically, we will consider a model where the 
dynamics of the Universe is driven in principle by pressureless matter and by the Chaplygin gas. The 
luminosity distance for the configuration where, in general, Chaplygin gas and dark matter are present, 
is evaluated from which the relation between the magnitude and the redshift z is established. The 
observational data are then considered, and they are fitted using four free parameters: the ratio of the 
density fraction, with respect to the critical density p c , of the pressureless matter and of the Chaplygin 
gas, fl m o and & c q respectively, the sound velocity of the Chaplygin gas A, in terms of the velocity of 
light, the curvature parameter £lko and the Hubble parameter Hq. All these parameters are evaluated 
today. 

The fact that we consider a model containing pressureless matter and the Chaplygin gas, implies that 
in principle we ignore the possibility that the Chaplygin gas could unify the description of dark energy and 
dark matter |11II12| . This unification is suggested by the fact that, at perturbative level, the Chaplygin 
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gas does not agglomerate at very large scales while it may agglomerate at small scales. The possible 
unification of both dark components of the matter content of the Universe through the Chaplygin gas 
has increased the interest on this new possible exotic fluid. A recent analysis of type la supernovae data 
correlated with the CGM ^B] used extensively this idea, employing the generalized equation of state J2J) 
in a flat Universe and excluding, ab initio, the possibility to have a dark matter component. In this work, 
the authors tried to restrict the possible values of the parameter a, but their results indicate a yet large 
range of possible values for this parameter. In the Ref. the authors studied, using type la supernovae 
data, essentially two scenarios: a flat Universe with the generalized Chaplygin gas and dark matter; a 
non-flat Universe with the ordinary Chaplygin gas and baryonic matter (with a fixed Q(, = 0.04). For 
the last case, the most favoured configuration was obtained for a closed Universe with tt c o = 1.27 and 
A = 0.78, using a \ 2 statistic. As it will be seen later, this result is consistent with the more general 
analysis we will perform here. 

The unification program using the Chaplygin gas has been recently criticized [201 121| . In our opinion 
this question remains an open issue, for the following reasons. In Ref. [3U|, the authors consider a 
model containing only the Chaplygin gas, and find a spectrum for the clustering of matter that is not 
in agreement with the 2dFGRS observations, unless the parameter a in J3J) is very small, a < 10~ a . 
However, we must not forget that there are baryons in the Universe. Even the quantity of baryons is 
small, it presences claims at least for a two fluid model, whose behaviour is in general quite different of 
a single fluid model [22] ■ Moreover, we may rise the question of what is really observed in the matter 
power spectrum. In Ref. the authors are more stringent by stating that, when the large scale 

structure data are crossed with the CMB data, the likelihood in the configuration space leads to a range 
of value for the a parameter that is very near zero and A ~ 1, i.e., essentially a cosmological constant 
model. However, the CMB analysis suffers from a high degeneracy in the space of the parameters, and 
the crossing of these data with those coming from type la supernovae must bring again the Chaplygin 
gas to scene. Here, we open the possibility to have dark matter present, not excluding at the same time 
the possibility to have a pure CGM. 

Our goal here, in contrast with the Refs. |181 119| . is to perform an extensive and comprehensive 
analysis of the problem of fitting the type la supernovae data using the ordinary Chaplygin gas given by 
Q. This restriction on the parameter a keeps us in contact with a string motivation for the Chaplygin 
gas. At the same time, our intention is to leave space for the presence of other fluids besides the Chaplygin 
gas as well as a spatial curvature term, constraining these parameters (together with the Chaplygin gas 
sound velocity and the Hubble parameter) only through the confrontation with the supernovae type la 
data. In this sense, the present work is more general than the previous ones, excepting for the restriction 
on the equation of state, restriction dictated by the theoretical motivation for the Chaplygin gas. A 
detailed description of the Bayesian statistics method of analysis is presented, as well as the meaning 
of the observational limits on the different free parameters of the model obtained through this method. 
A comparison with the cosmological constant model (AC DM) is exhibited, which is one of the main 
purpose of this work. 

It will be shown, in particular, that the CGM , in what concerns the supernovae type la data, favours 
strongly, at la, a closed Universe (flko = — 0-84^]^ '| 3 ), peaked in the zero value for the cold dark matter 
density (Q m o = O.O^q 2 ) and a sound velocity c 2 s = 0.93+q;2iC near, but not equal, to the value that 
corresponds to the cosmological constant. The present age of the Universe is to = 14.2^'g Gy in the 
CGM, a bit smaller than that in the ACDM model (to = 15.4^'gGy), but still compatible with other 
astronomical measurements The deceleration parameter is highly negative: qo = — 0.98+J'g2- The 
CGM is always prefered with respect to the ACDM model, with a confidence level of 55.3% when 
all free parameters are considered. This confidence level is considerably higher if the curvature or the 
dark matter density is fixed to zero. An almost flat Universe is predicted if dark mater is absent, 
^fco = 0.17j^'5g(f2 m o = 0). Later, we will make a brief comparison between the results obtained with 
the recent data coming from the anisotropy of the cosmic microwave background (including the data 
coming from the WMAP) |24l I25j . 2dFGRS measurements j2B| and Lyman a forest even if a proper 
combined analysis is outside the scope of the present paper. 

In next section, the Chaplygin gas model is described. In section 3, the Bayesian probabilistic analysis 
is described. In section 4, the cosmological parameters are estimated, and in section 5 our conclusions 
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are presented. 



2 The Chaplygin gas model and the least-squares fitting 

We now proceed by constructing the cosmological model based on the Chaplygin gas, which we will test 
against the type la supernovae observations. We must specify the observable cosmological parameters in 
the model: we will deal with the Hubble parameter Hq, the curvature parameter £lko, the dark matter 
density parameter f2 m o, the Chaplygin gas density parameter f2 c o and the Chaplygin gas sound velocity 
Cg. All these parameters are evaluate today. 

The sound velocity of the Chaplygin gas today, in units of c, is given by 

cl=A = 4- ■ (4) 

PcO 

The ratios of the density fractions with respect to the current critical density are related by ilko + fi TO + 
fi c = 1. For the Chaplygin gas case we take O m o ^ and Sl c o ^ 0, so flko ^ 1- When the limit A — > 1 is 
taken (c = 1 hereafter), the cosmological constant model is obtained from the Chaplygin model, so f2 c o is 
relabeled as Q\ (the ratio of the vacuum energy density with respect to the current critical density), and 
the only range restriction is £l m o ^ 0. It will be verified that the best-fitting and the Bayesian inference 
suggest flmo ~ 0, this result becomes quite interesting if we take into account some recent considerations 
about a unification of cold dark matter and dark energy in Chaplygin gas models [T5|. Moreover, closed 
models will be favoured. 

The equation governing the evolution of our model is 

/ D~\ 

(5) 

\Ub I LL~ O \ U~ V IX" / 

It can be rewritten as 

/a\ 2 I Cl . I i T o._\ 

(6) 

/ 

where Hq is the Hubble parameter today f2j = Bi. (^ denoting a matter component and p c the critical 

density), ilko — §^|j and the scale factor was normalized to unity today a = 1. 

The luminosity distance is obtained by standard procedures |28| . using the equation for the light 
trajectory in the above specified background, and its definition, 

1/2 

4tt I 




^=lif) (7) 



where L is the absolute luminosity of the source, and I is the luminosity measured by the observer. This 
expression can be rewritten as 

D L = (l + z)r , (8) 

r being the co-moving distance of the source. Taking into account the definitions of absolute and apparent 
magnitudes in terms of the luminosity L and I, M and m respectively, we finally obtain the relation valid 
for the three types of spatial section: 

c 1 



M = 51ogUl + z)— -==></ V\Sho\x (9) 



x 



dz' 



^fWi + z') 3 + n cX) ^A + 0~- A)(i + z'f + n fc0 (i + z'f 



(10) 



where f(x) = sinh(cc) for Sl^o > (open Universe with k < 0), f(x) — sin(x) for Clko < (closed Universe 
with k > 0) and f(x) — x for O^q = (flat Universe with k = 0). 
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Following the same lines, we can obtain the expression of the age of the Universe today: 



r t . (id 

which will be another parameter evaluated for the different models. 

The deceleration parameter q is defined as q = — 1 — H/H 2 , and its value calculated today, using Eq. 
®), reads 

90 = ^ , (12) 

so depending on the values of the three parameters above, an accelerating Universe (qo < 0) or decelerating 
Universe (go > 0) is obtained. When A = 1 the well-known result for the cosmological constant model is 
recovered, q a = (Cl m0 /2) - fl c0 . 

Another useful result is the calculation whether the Universe will expand eternally or not from the 
present time. A simple method used here searches for the roots of a(t) in the r.h.s. of Eq. ©; if there is 
one that is a real value and greater than ao = 1, then the Universe is not always expanding in the future. 

We proceed by fitting the SNe la data using the model described above. Essentially, we compute the 
quantity distance moduli, 

Mo = 5 log (^] + 25 , (13) 



\Mpc / 

and compare the same distance moduli as obtained from observations. The quality of the fitting is 
characterized by the % 2 parameter of the least-squares statistic, as defined in Ref. [2], 

i fiQ,i tnz,i 

In this expression, /Ltg i is the measured value, i is the value calculated through the model described 
above, er^ o i is the measurement error, a^ nz i is the dispersion in the distance modulus due to the dispersion 
in galaxy redshift due to peculiar velocities. This quantity we will taken as 

d\ogD L 

O-mz = ^ Oz , (15) 

oz 

where, following Ref. 0EH], a z = 200 km/ s. 

In table H the data concerning 26 type la supernovae with the error bar are displayed. These are 
essentially the type la supernovae employed in the first works on the problem of the acceleration of 
the Universe 0. We restrict ourselves to this sample (samples with up to about one hundred type la 
supernovae are now available), since it contains some of the better studied SNe la. It must be stressed 
that since one of the goals of this work is to compare competitive models, the choice of the sample is not 
so essential, provided it is not too small neither contains doubtful data. 

There are two methods to determine the relationship between the shape of SNe la light curve and its 
peak luminosity: MLCS (Multicolor Light Curve Shapes) [35] and a template fitting method (Atoi5(£?)) 
|33|. For both methods, the fit determines the light curve parameters and their uncertainties. The set of 
26 SNe la used in this article are within the expected statistical uncertainties range of the two methods. 
Specifically, we use the parameters obtained by the template fitting method in our analyses. 

In the tables [21 and we evaluate, in fact, \l, the estimated errors for degree of freedom, i.e., % 2 
divided by 26, the number of type la supernovae chosen in this article. The values of xt f° r the best- 
fitting are indeed small, between 0.75 and 0.78, if compared to other SNe la analyses 0, due to our choice 
of SNe la that avoided those with large observational errors. In previous works jHOll^, the restricted 
case of spatially flat Universe was analysed by means of best-fitting parameter estimation. 

The table lists the best-fitting parameters of the cosmological constant model (the limit A — ► 1 
of the Chaplygin gas model). When the spatial section is open or flat, the best- fitting favours a flat 
Universe with Hq = 61.8 km/Mpc.s, f2 m o/^A = 0.24/0.76 and to = 16.5 Gy, approximately the same 
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0.014 


34.13 (0.14) 


1992bp 


0.080 


37.96 (0.15) 
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0.018 


34.88 (0.21) 


1992br 


0.087 


38.09 (0.36) 
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0.020 


34.77 (0.15) 


1992aq 


0.111 


38.33 (0.23) 


1992ag 


0.026 


35.53 (0.20) 


1996J 


0.30 


40.99 (0.25) 


1992P 


0.026 


35.59 (0.16) 


1996K 


0.38 


42.21 (0.18) 


1992bg 


0.035 


36.49 (0.21) 


1996E 


0.43 


42.03 (0.22) 


1992bl 


0.043 


36.53 (0.20) 


1996U 


0.43 


42.34 (0.17) 


1992bh 


0.045 


36.87 (0.17) 
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0.44 


42.26 (0.16) 
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1993ag 


0.050 


37.11 (0.19) 


1997cj 


0.50 


42.70 (0.16) 


19930 


0.052 


37.31 (0.14) 


19961 


0.57 


42.83 (0.21) 


1992bs 


0.064 


37.63 (0.18) 


1996H 


0.62 


43.01 (0.15) 


1992ae 


0.075 


37.77 (0.19) 


1997ck 


0.97 


44.30 (0.19) 



Table 1: The SNe la data of the 26 supernovae used in this article, obtained by the template fitting 
method (Amis (£?)). 
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model with k < 


Cosmological constant 
model with k > 




0.7743 


0.7539 
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61.8 


62.4 


^fcO 


0.00 


-0.80 




0.24 


0.57 




0.76 


1.23 


to 


16.5 


15.7 


qo 


-0.64 


-0.95 



Table 2: The best-fitting parameters, i.e., when xl 1S minimum, for each type of space section of the 
cosmological constant model. Ho is given in km/Mpc.s and to in Gy. 

result estimated by other research groups [3 EH- But a closed Universe gives a better fitting of the 
parameters, since xt is smaller, suggesting a different Universe, with positive curvature and younger. It 
is clear that small amounts of xt lead to totally different parameters, i.e., the parameter estimation is 
highly dependent on xt- 

The best- fitting parameters for the Chaplygin gas model are given by table |21 Slightly smaller values 
of xt favour the Chaplygin gas model over the cosmological constant model, independent of the spatial 
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Parameters 


Chaplygin gas model 
with k < 


Chaplygin gas model 
with k > 
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0.7606 


0.7527 


Ho 


62.1 


62.5 




0.00 


-0.44 
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n nn 
u.uu 


U. lo 


FlcO 


1.00 


1.31 


A 


0.847 


0.858 


to 


15.1 


14.6 


qo 


-0.77 


-0.97 



Table 3: The best-fitting parameters, i.e., when xl is minimum, for each type of space section of the 
Chaplygin gas model. H is given in km/Mpc.s, A in units of c and t in Gy. 

section type. The best-of-all fitting suggests a closed Universe dominated by the Chaplygin gas, featuring 
a positive curvature, small value of Q m o, and the smallest age. In the case of open or flat spatial section, 
the best-fit gives a limit case of a flat Universe without cold dark matter content, just filled by the energy 
density of the Chaplygin gas. 

Nevertheless, the xl best-fitting analysis has many limitations; the following questions cannot be 
answered. How much worse is one parameter set compared to other one, for example in table [21 or [21 ? 
What is the likelihood of a closed Universe using the Chaplygin gas model (table |3J) ? And the likelihood 
of the age of Universe to be in the range 13 to 15 Gy 1 

3 Bayesian probabilities 

We can answer these questions using the frequentist (traditional) or the Bayesian probability theory of 
statistics. See Refs. |^[3Sl[3niEZl f° r discussions about them and some applications in physics. We have 
chosen the Bayesian approach to avoid some problems of the frequentist theory: the need of simulating 
the observational data, lack of mathematical formalism, etc. 

The posterior probability density function (PDF) for the cosmological parameters Hq, fi m o, Q c o and A, 
given the set of distance moduli fio data and a cosmological model, can be obtained from the probability 
of the no data conditional on (H , fi m o, ^ c o> A), if we use the Bayes's theorem: 

iv o r> j \ KMo I #o,fi m0 ,Q c o,A) p(H ,£lmo,n c0 ,A) 

p[Ho,\l m n,\l c0 ,A \ Hq) = — , (loj 

PifJ-o) 

where the prior probabilities p(Ho, fi m o, ^co, A) and p(po) are considered here to be constants because 
we have no prior constraints on the parameters or on the data, besides some forbidden regions of the pa- 
rameter space corresponding to unphysical universes. So, the probability of (-ffo? ^mOi ^co, A) conditional 
on fio is proportional to the probability of /jq conditional on (Hq, f2 m o ; ^c0: A): 

p{H , Cl m0 , tt c0 , A\[io)<x p(po I H Q , Cl m0 , n c0 , A). (17) 

Each distance moduli fio,i is considered as independent with a Gaussian distribution, then the prob- 
ability of the set of distance moduli fio conditional on (Ho, n m o, f2 c o, A) is the product of the Gaussians, 
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as explained in Ref. 0: 



p((j, | H , Cl m0 , tt c0 , A) = Y[ I exp 

* V 27r Ko,i + (T L,i) 



(18) 



which can be written using the x 2 calculated previously, Eq. I|14fl . because the product of exponentials 
is the exponential of the sum present in x 2 , 

p(fi | H ,n m0 ,n c0 ,A) oc exp (^y) ' ^ 19 ^ 

Finally, by using the relation i|17|) . the probability of the parameters (Ho, f2 m o, Q c o, A) conditional on 
the set of distance modulus /io can be written in a normalized form: 

exp(-4H 

p(H ,n m0 ,n c0 ,A | no) = i L (20) 

Too di?o JZo d^™o d^co r« «p (- £ ) 

where the integrals are performed on the allowed region of the parameter space. 

Obviously, as the least-squares error x 2 statistic is minimum, the probability is maximum. The values 
of the parameters (H Q , £l m o, ^ c o, A) of this maximum are the most likely parameter values, when thought 
simultaneously. 

But when we ask for the best value of one parameter, the \ 2 analysis is not robust. For example, 
the x 2 minimum can be located in a narrow region with small values of x 2 whereas another larger region 
can also have small values of x 2 > an d this information is not taken into account to determine one of the 
parameters (Hq, Q m o, Qco, A) for the x 2 minimum. 

With the Bayesian probability i|2U|l . or likelihood, or PDF, defined in the four-dimensional parameter 
space, we can construct a probability for one parameter by marginalizing, i.e., integrating with relation 
to the other parameters, the so-called marginal probability. For example, the marginal probability for 
H is defined as 



p(H Q \no)= dQ m0 dfl c0 p(H ,n m0 ,n c o,A\ fi a )dA, (21) 

J —oo J —oc J — oo 

where again the integration region is restricted to the parameter space with physical meaning. The peak 
of this PDF provides the maximum likelihood estimate of the parameter Ho, which is usually different 
from the value of Hq for x 2 minimum, because generally, p(Hq, Q m o, £l c o, A | /zo) is not a four-dimensional 
Gaussian PDF (or some distribution alike). This maximum likelihood estimate of Hq is a more robust 
estimate because the rest of the parameter space is the integrate region for each value of Hq, so a larger 
region with high PDF values (but not the maximum) is not discarded. 

The availability of a PDF for a parameter also allows us to calculate the likelihood of arbitrary 
hypothesis, for example: H greater than some value, H between some range, etc. They are performed 
by calculating the CDF (cumulative distribution function), i.e., the integral of the PDF over the specified 
region. The inverse problem also appears: for a given likelihood value, calculate the region in the 
parameter space which has a CDF equal to some likelihood value. 

A well-known example is the credible or confidence region of la, 2a and 3tr likelihoods. To define 
this type of region, we use the property that Gaussian PDF function, cr _1 (27r) -1 / 2 exp(— x 2 /2a 2 ), when 
integrated over the region — la < x < la, gives a probability of approx. 68.27%; with the region 
—2a < x < 2a the obtained likelihood is approx. 95.45%; and the range —3a < x < 3a yields a CDF of 
approx. 99.73%. So, a values mean probabilities. 

For an arbitrary PDF function p(x), the calculation of the 2a (for example) credible or confidence 
region in one dimension is usually defined as follows: obtain the PDF level between zero and the peak 
PDF (when x = xq) which have the intersection with p(x), X- and x + , so that f*_p(x)dx ~ 95.45. The 

estimation of x using 2a credible region is described as (^o)xt-^oJ meaning that the PDF is peaked at 
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x and the CDF of the region x_ < x < x + is equal to 2a (95.45%). See, for example, the thin line 
PDF in figure^] with its la and 2a credible regions, PDF levels and intersections. In n dimensions, the 
intersection produces a n-dimensional region that becomes the integrate region, for an example in two 
dimensions, figure 01 shows the la, 2a and 3a credible regions. 

The marginal joint probability (PDF) as function of two parameters, for example p(£l m o, ^co), is an 
integral in the two dimensional parameter space of (Ho, A), 

/oo />oo 
Hq / p(Ho,il m Q,n c0 ,A\ fi )dA. (22) 
-oo J —oo 

Analogously, the marginalization method can be used to obtain, for example, p{Hq, fi m o, H c o | no) from 
p{Hq, f2 m o, O c o, A | /iq) by means of one integral over the A parameter space (or marginalizing over A). 

Any quantity depending on the parameters can also have a probability for it. For example, in the 
case of to, the dynamical age of Universe, the PDF is obtained from: 

/OO />00 />OG pOO 

dH / d£! m0 / d£! c0 / p(H , ^mo, ^co, A | t , Mo)<b4, (23) 
-oo J —oo J —oo J —oo 

where the integrand gives the likelihood of the parameters which give a certain age to, so the four- 
dimensional integral is a sum of the likelihood of all possible parameters which produces a Universe with 
age to ■ To avoid the computation of a four-dimensional integral for each value to , it is better to calculate 
the age of Universe today for the four-dimensional parameter space and store the cumulative probabilities 
for each value of to ■ 



4 Analyses of the estimated parameters 

We have performed some long calculations using the Bayesian approach to obtain the parameter estima- 
tions and answers for some hypothesis. In the Chaplygin gas case, the five parameters, {Hq, &ko, ^mOj ^co, A), 
the age of the Universe to and the deceleration parameter qo were estimated with a central value and 
a 2a (95.45%) credible region. Each independent parameter estimate used a marginal likelihood of the 
type of Eq. 1)21(1 . where three-dimensional integrals are computed for each value of the parameter (and 
in the integrand, % 2 needs the calculation of about one hundred numerical integrals, or four times the 
number of supernovae). 

An ideal calculation to compute the n-dimensional integrals would include an infinite number of 
samples of a parameter space with infinite volume, but in practical estimations we chose a finite region 
of the parameter space (such that outside it the probabilities are almost null) and a finite number of 
samples. Considering n samples for each parameter dimension, the estimation of one parameter needs 
n 4 computations of x 2 ( as there are four independent parameters), and by proper marginalization the 
PDF of each parameter is calculated. For this article, n « 30, and the number of computations of \ 2 
is at least n A ks 10 6 just for the Chaplygin gas model (some calculations are repeated because it is not 
worth storing 10 6 results due to computational memory constraints), or 100n 4 sa 10 8 numerical integrals 
of the type in Eq. (|10|1 . We also calculate n 4 « 10 6 times the age of Universe today to discard paramater 
space regions which are not physial (to is not real, etc) and to estimate to- Likewise, the qo estimate 
demands n 4 calculations of Eq. i|12|) f° r Qo- The probability of eternally expanding Universe consists 
of n 4 calculations of the roots of Eq. © . The other cases listed in tables 0] and [S] demands n 3 or n 2 
calculations, instead of n 4 . 

Tables 0] and summarizes the Bayesian parameter estimation results for seven different cases. As 
expected, the central values of the Bayesian parameter estimation and the best-fitting parameters of 
tables |2 or |3 are different. In most cases, the likelihood functions (PDF) of the parameters have a shape 
similar to a Gaussian function, but with some asymmetry, consequently they are not shown here because 
the central value and credible region description are enough to realize the PDF behaviour. 
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I -) 1 

Parameters 


Cosmological constant 
model 


Cosmological constant 
model with fc = 


Cosmological constant 
model with i2 m0 = 




62.2 ± 3.1 


62.2 ± 3.1 


<ii /( 1 OA 

61.4 ± 3.0 


o 

"feO 


-0.80_ 1 34 


A 

o 


A cc+0.50 

0.55_ 037 




n co+0.56 
°- 58 -0.58 


A r, /I +0.21 

0- 24 -o.i6 







1 r>1 +0 81 

!- 21 -0.91 


n -^+0 16 

0.76± .2i 


„ . r -i_o 37 

0-45-0.50 


to 


1 r ,1+3 4 


15.51 3 ; 3 


19.3_3 i 


9o 


n QQ+0-75 
u - ay -0.52 


_U - '0-0.14 


n 47 +0.50 
_U - 4 '-0.37 


p(fl k0 < 0) 


81.6%(1.33cr) 




0.03% 


p(^A > 0) 


99.8%(3.03cr) 


99.999993%(4.86cr) 


95.7%(2.02cr) 


p(qa < o) 


99.6%(2.86cr) 


99.97%(3.64cr) 


96.1%(2.06ct) 


p(ct > 0) 


93.7%(1.86cr) 


99.1%(2.61cr) 


37.8%(0.49cr) 



Table 4: The estimated parameters for the cosmological constant model, using the Bayesian analysis to 
obtain the peak of the marginal probability and the credible region for each parameter. H is given in 
km/Mpc.s, A in units of c and to in Gy. 



Parameters 


Chaplygin gas 
model 


CGM 
with k = 


CGM 
with Q m o = 


CGM with k = 
and f2 TO o = 0.04 


Ho 


62.1+ji 


61.4 ±2.8 


61.9 ±2.8 


61.8 ±2.8 


^fcO 







17+ - 83 

' -1.58 







0-00i°;^ 


u - uu -o.oo 





0.04 


n c0 > o 


1 40+ 1 - 15 


i n +o.oo 
^"-o.ss 




0.96 


A sc 1 


Q3+ 07 


93+ 07 


0-78±g;?l 


o.87ig : g 


to 


14.2±?;1 


i3.i±?;S 


14.51?;? 


14.8l?-j 


qo 


-0.98l ;°J 


-0.65l°;f 7 


C S6+ - 61 
U.OD_ 26 


-o.8i±8-.g 


p{^w < o) 


84.0%(1.41cr) 




55.7%(0.77cr) 




p(qo < o) 


97.5%(2.24ct) 


99.96%(3.55cr) 


95.2%(1.98cr) 


99.997%(4.17cr) 


p(a > 0) 


93.9%(1.87<t) 


99.7%(2.99cr) 


63.8%(0.91cr) 


100% 



Table 5: The estimated parameters for the Chaplygin gas model, using the Bayesian analysis to obtain 
the peak of the marginal probability and the credible region for each parameter. Ho is given in km/Mpc.s, 
A in units of c and to in Gy. 



10 



4.1 The Hubble parameter H 

The parameter estimation of Hq is almost the same for all cases, around 62 km/Mpc.s with a narrow 
credible region of ±3 km/Mpc.s. This estimation is compatible with other SNe la analyses [21 E]. Hence, 
the choice of Chaplygin gas model or KCDM does not have any significant consequence on the Hq 
estimates we have made. 

4.2 The curvature density parameter Qko 

Likewise, the curvature density parameter O^o estimation is almost the same, — 0.84l*;§3 and — 0.80^' H, 
assuming the Chaplygin gas model or KCDM, respectively, and in both cases a closed Universe (fifco < 0) 
is preferred at 84.0% (1.41<r) and 81.6% (1.33cr) confidence levels, also respectively, in agreement with 
the conclusions of many other previous works on the subject OEj. In the CGM , a spatially flat Universe 
is ruled out at 67.5% (0.98a) confidence level, i.e., the PDF of fi^o = is smaller than the PDF of 
-1.42 < Q k0 < and this region has a CDF of 67.5% (0.98er). The KCDM model shows a similar 
behaviour, a spatially flat Universe is ruled out at 63.3% (0.90cr) level, i.e., the PDF of flko — is smaller 
than the PDF of -1.46 < fi fc0 < 0. 

But assuming the hypothesis that fi TO o = 0, describing a Universe empty of cold dark matter and 
totally dominated by the Chaplygin gas, we have quite different estimates. For the CGM, the Bayesian 
analysis gives flko = O-^tiJg, with a closed Universe favoured at 55.7% (0.77<r), and a flat Universe is 
strongly preferred at 83.3% (1.38cx) confidence level, i.e., the PDF of flto = is smaller than the PDF of 
< flko < 0.30 and this region has a small CDF of 16.7%. 

The same hypothesis, fl m o — 0, applied to the KCDM model, leads to the estimation Q^o = 0.55^q'Ij, 
so now an open Universe is quite favourable at 99.97% (3.60cr) confidence level. A spatially flat Universe 
is ruled out at 99.92% (3.36er) confidence level, i.e., the PDF of flko = is smaller than the PDF of 

< n k0 < 1.43. 

Therefore, the O m o = case clearly discriminates the Chaplygin gas and the KCDM models, as their 
estimates are quite different. This discrimination is still enhanced if the additional hypothesis of spatially 
flat Universe is assumed (to be compatible with inflationary models of the primordial Universe, or some 
CMB estimations, for example), when the CGM agrees at a high confidence level while KCDM is ruled 
out. 

4.3 The cold dark matter density parameter f2 m0 

But the cold dark matter density parameter fl m o depends quite on the cosmological models. Figure ^ 
clearly shows that the PDF for the Chaplygin gas and the KCDM models behave differently, fl m o = 
O.OOto'oo an d ^mo = 0.58lg'gg, respectively. So the Chaplygin gas model favours small values of f2 m0 
at a high level of confidence, for example, fl m o < 0.39 at la (68.27%) confidence level, while KCDM 
estimates fl m o < 0.39 with a probability of 31.1% (0.40cr), or, Q m0 = 0. 58^28 a * ^ a level. 

While the most favoured value of fl m o is zero for the CGM, the cosmological constant model rules 
out f2 m o = at 89.5% (1.62er) confidence level, i.e., the PDF of £l m o = (equal to 0.56) is smaller than 
the PDF of < n m0 < 0.99 which has a CDF of 89.5%. 

On the hypothesis that the Universe is spatially flat, the CGM and KCDM have narrower confidence 
regions, see tables |2 and 03 It is worth noting that the CGM continues to indicate a fi m0 peaked at 
zero, but the KCDM model rules out fl m o = completly, at 99.94% (3.44cr) confidence level (when 
< O m0 < 0.63 the PDF is greater than the PDF of Q m0 = 0). 

Based on the SNe la data, we can conclude that, assuming the Chaplygin gas model, the estimated 
values of f2 m o are decreased (with respect to KCDM), or if independent estimations of il m o suggest low 
values then the Chaplygin gas model is favoured over the KCDM model. 

4.4 The Chaplygin gas density parameter £7 c0 

The Chaplygin gas density parameter f2 c o estimation is quite spread, i.e., it has a large credible region, 
fl c o = 1.40]ti'i6- O nrv f° r t ne particular case of a spatially flat Universe (k = 0) the estimation of fi c o is 
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Figure 1: The PDF of Q m o for the Chaplygin gas model (bold lines) and the cosmological constant model (thin lines). 
The solid lines arc the PDF, the 2cr (95.45%) credible regions are given by dashed lines and the lc (68.27%) regions are 
delimited by dotted lines. The Chaplygin gas case has maximum probability 2.55 for Q m o = 0.0, the 2<r credible region 
reads $C QmO < 0.82 with probability level 0.26, and the lc credible region is given by ^ H m o < 0.39 with PDF level 
1.10. The ACDM has probability peak of 1.11 when Q m o = 0.58, the 2<r credible region reads ^ f2 m o < 113 with 
probability level 0.31, and the lcr credible region is given by 0.20 $C f2 m () < 0.88 with PDF level 0.80. This graphics clearly 
shows that small values of Q m o are more preferred by the Chaplygin gas model. 

less spread. The f2 m o = case gives a lower estimate but still strongly spread: fl c o = 0.83^'^. 

For the ACDM model obtained from the Chaplygin gas model as A = 1, the dark energy density f^A 
is positive at a 3.3c (99.8%) confidence level, and the estimation Ha = l-2llo'gi becomes more peaked 
for a spatially flat Universe, Ha = 0.76^g'2i- See table[5]for more estimates. 

4.5 A, the sound velocity of the Chaplygin gas 

The analysis of the parameter A is very conclusive, because ACDM is a special case of the Chaplygin gas 
model when A = 1, therefore we can compare the probabilities of each cosmological model. Indeed, figure 
[5]shows that the best estimation reads A = 0.93^0 21: with peak probability level 4.30, while the ACDM 
limit (A = 1) has a slightly lower likelihood of 4.12. More specifically, the Chaplygin gas model is the 
more preferred model over the range 0.87 ^ A < 1 because the PDF is greater than 4.12, and this region 
represents 55.3% (0.76er) of the total probability. So, even if there is no additional hypothesis about the 
curvature parameter £lko or the cold dark matter density parameter O m o, the Bayesian analysis explicitly 
estimates that the Chaplygin gas model is more likely than ACDM at 0.76cr (55.3%) confidence level, 
although the peak likelihood is just 4% greater than ACDM likelihood. 

Under the hypothesis that £l m o = 0, so the Chaplygin gas is the only content of the Universe, figure 
El gives A = 0.781q29 an d indicates that the Chaplygin gas model has a peak likelihood of 4.86, quite 
greater than the ACDM model probability level of 0.87 (for A = 1), and on a broader region 0.67 ^ A < 1 
the PDF is greater than 0.87. The parameter A is inside this region with likelihood of 91.8% (1.74cr), 
i.e., the Chaplygin gas model is preferred over the ACDM model at 1.74cr (91.8%) confidence level. 

Assuming now the hypothesis of a flat Universe, Sl^o = 0, we obtain A = 0.93t^"o 

with a peak 
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Figure 2: The PDF of A for the Chaplygin gas model (bold lines) and the same with H m o = (thin lines). The solid 
lines are the PDF, the 2a (95.45%) credible regions are given by dashed lines and the la (68.27%) regions are delimited by 
dotted lines. The non-restricted case has maximum probability 4.30 for A = 0.93, the 2a credible region reads 0.72 ^ A < 1 
with probability level 0.85, and the la credible region is given by 0.84 $C A < 1 with PDF level 3.77. The probability 
of the ACDM obtained when A = 1 is 4.12, so the Chaplygin gas model is more likely (has greater PDF) in the range 
0.87 ^ A < 1, which has a CDF of 55.3% (0.76cr). The f2 m0 = case has probability peak of 4.86 when A = 0.78, the 2a 
credible region reads 0.59 SC A < 1 with probability level 0.25, the la credible region is given by 0.71 SC A < 0.89 with PDF 
level 2.39, and the ACDM has a PDF level of only 0.87 which suggests that the Chaplygin gas model is strongly preferred 
if there is no cold dark matter (f2 m o = 0) content in the Universe, because the range 0.67 SC A < 1 has a CDF of 91.8% 
(1.74<t) such that the PDF level is greater than 0.87. 

likelihood of 5.12, greater than the ACDM model probability level of 4.18 (for A = 1). The region 
0.87 s: A < 1 has a PDF greater than 4.18, which has a CDF of 63.5% (0.91a-), i.e., the Chaplygin gas 
model is preferred over ACDM at 63.5% (0.91<r) confidence level. 

In the framework of the cosmological model estimations based on SNe la, the baryonic matter has the 
same behaviour of the cold dark matter, so we can assume a typical value of 0.04 for the baryonic density 
parameter by setting fl m o — 0.04. Still imposing the hypothesis of a flat Universe, the last column of 
table |3 shows that A — 0.87^q']J with a maximum PDF of 5.56 and the region 0.59 ^ A < 1 having 
greater PDF than the PDF of ACDM limit (A = 1), 0.13. The CDF of this region means that the 
CGM is preferred over ACDM at 99.4% (2.72c) confidence level, under the hypothesis £lko — and 
!l m o = 0.04. Not shown in [2 the ACDM model parameter estimates read, Ho — 64.4^2 \ km/Mpc.s 
and t = 24.2±° ;| Gy, with an age of the Universe excessively higher than those obtained through other 
astronomical estimations. 

Combining both hypothesis, £lko = and f2 m o = 0, the parameter estimation is quite physically 
acceptable for the CGM. Not show explicitly in table we have obtained Hq = 61.9^2 ff km/Mpc.s, 
A = 0.85±g;5J, t = U.7tf 7 5 Gy and q = -0.81±g;fl with p(q a < 0) = 99.996%(4.12cr). The Chaplygin 
gas model is strongly favoured over the ACDM model at 99.93% (3.38er) confidence level, as the PDF 
of A — 0.85 is equal to 5.53 and the likelihood of A = 1 is 0.02, so a large region 0.47 ^ A < 1 
has PDF levels greater than 0.02. Moreover, the ACDM model produces totally unphysical estimates: 
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to = 3681^ x 10 3 Gy. 

4.6 The present age t of the Universe 

Looking at the dynamical age of the Universe today, to , given in tables 0] and we verify that the 
Chaplygin gas model estimates lower values of to as well as narrower credible regions (to — 14.2j^'g), 
when compared to the estimations of the ACDM model (to — 15.4^'g). By means of independent 
estimations of to from different observations (globular agglomerates, etc), the Chaplygin gas model would 
be more preferred if, for example, these estimations give to between 13Gy to 14Gy. 

If we impose some constraints, k = or O m o = 0, we still have lower estimates for to of the CGM 
than the ACDM model estimations. The latter are strongly sensible to the value of £l m o'- for example, 
if k = 0, we obtain t Q = 24.2 ±° ;| Gy for fl mQ = 0.04 and totally unphysical value of t = 368±^ x 10 3 Gy 
for Q m o = 0. Quite different, the CGM produces to values with a weak dependence on Qko and £l m a- 

4.7 The deceleration parameter q of the Universe 

All the estimations based on the Chaplygin gas model or the cosmological constant model suggest an 
accelerating Universe today, i.e., <?o < 0, with a high confidence level (at least more than 95%). 

4.8 The eternally expanding Universe 

The probability of having a eternally expanding Universe, given here by demanding a > during the 
future evolution of the Universe, is estimated to be almost one (at least more than 93%), excepting in 
the fl m o — case for both the CGM and ACDM when p(a > 0) is far from unity, see tables 0] and [5] 

4.9 Two-dimensional analysis in the (Q m o,Q A ) and (Q m o,Q c o) parameter space 

Finally, it is useful to analyze the two-dimensional joint PDF, for example, in the parameter space of 
(Q m o,£l A ), where the PDF p(r2 m0 ,^A | Mo) is given by the integral of p(H ,^ m o,^A I Mo) over the H 
space. The well-known credible regions for the ACDM model are reproduced here, see figure [3J It is 
worth noting that the probability peak gives (f2 m 0j Oa) = (0.60, 1.27) which are different from the central 
values listed in table £l m o = 0.58 and ft a — 1.21, due to the integration over the 17a space to obtain 
p(Q m o | fio) from p(Q m Q, J7a | Mo): and analogously for p(fl\ | Mo)- The interpretation is the following : 
(fi m o, = (0.60, 1.27) means that these values are the most likely simultaneously values of (Cl m o, ^a) 
independent of the parameter Hq] while Q m o = 0.58 is the most likely value of f2 m o independent of the 
other parameters (Ho,Qa), and fl\ = 1.21 is the most likely value independent of the other parameters 
(Ho,Q m o). Figure ^ shows p(Q rn o \ (1q) which is obtained by integrating p(Q rn o, J7a | Mo) over the £l\ 
space. 

The contour curves are defined such that the enclosed regions have a cumulative probability equal 
to la (68,27%), 2a (95,45%) and 3a (99,73%) credible levels. For example, this means simply that 
fl m o and Qa have values simultaneously inside the 2a (95,45%) region with a likelihood of 95,45%. The 
line f2 m o + = 1 represents a spatially flat Universe (k = 0, £lko — 0), the parameter space above 
corresponds to a closed Universe (fc > 0, ttko < 0), and below, to an open Universe (fc < 0, Qko > 0). As 
the closed Universe region clearly dominates and presents high PDF values, a closed Universe is estimated 
at 1.33cr (81.6%) confidence level. 

For the case of the Chaplygin gas model, the PDF p(Q m o, fl c o Mo) is a two-dimensional integral given 
by Eq. (|22|l . and figure 01 displays its behaviour with credible regions quite different from the ACDM. 
The PDF peak is now located at (f2 m o, ^co) = (0.0, 0.85), contrasting with the central values fl m o = 0.00 
and il c o = 1.40 of table |5j due to the Bayesian integrations. Of course, the integral of p(f2 m o, f2 c o | Mo) 
over the f2 c o space yields p(£l m0 | /i ), shown in figure ^ Once more, a spatially closed Universe is 
strongly favoured at 1.41(7 (84.0%) confidence level, despite the PDF peak being inside the spatially open 
Universe region. 

These two graphics also emphasize the totally different likelihoods for (fifco,^mo) = (0,0.04) in the 
CGM and ACDM model. This point corresponds to (fl m0 ,n A ) = (0.04,0.96) in figure for ACDM, 
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Figure 3: The graphics of the joint PDF as function of (O m 0;^A) for the ACDM, where p(f2 m o,f2A) is a integral of 
p(Hq, O m 0i ^a) over the Hq parameter space. It shows the well-known shape of the confidence regions for the cosmological 
constant model. The joint PDF peak has the value 2.07 for (Q m o,Q/^) = (0.60,1.27) (shown by the large dot), the la 
(68,27%), 2a (95,45%) and 3a (99,73%) credible regions have joint PDF levels of 0.85, 0.15 and 0.01, respectively. As 
QlcO + ^mO + ^A = 1, the probability for a spatially flat Universe is on the line Q m o + ^A = 1, above it we have the region 
for a closed Universe (k > 0, fl^o < 0), and below, the region for an open Universe (k < 0, H^q > 0), so the integral of the 
PDF over the upper region gives the likelihood of a closed Universe being 81.6%. 

with a PDF level of 0.02, well below the PDF peak 2.07, so a large region of the (£l m o, ^a) parameter 
space has greater PDF than 0.02 yielding a CDF of 99.6% (2.91a), i.e., the simultaneously hypothesis 
of flat Universe and il m o = 0.04 (representing a typical baryonic density parameter of 0.04) is ruled out 
at 99.6% (2.91er) assuming the ACDM. The limit case of (ilko,Q m o) = (0,0) gives a confidence level of 
99.93% (3.39cr). 

Figure 0] clearly shows that the point (f2 m o,r2 c o) = (0.04,0.96) is near the maxima of PDF, their 
PDF levels are 1.41 and 1.47, respectively. Therefore, the region with smaller PDF than 1.41 is almost 
the whole of the parameter space, and its CDF is equal to 98.6%, i.e, for the CDM the simultaneously 
hypothesis of flat Universe and f2 m o = 0.04 is favoured at 98.6% (2.45cr) confidence level. In the same 
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Figure 4: The graphics of the joint PDF as function of (Q m o,n c o) for the Chaplygin gas model, where p(Q m Q,Q c o) is a 
integral of p(Hq, Q m o, O c q, A) over the (Ho, A) parameter space. It shows a different shape for the confidence regions with 
respect to the cosmological constant model. The joint PDF peak has the value 1.47 for (Q m o,Q c o) = (0.00,0.84) (shown 
by the large dot), the lcr (68,27%), 2a (95,45%) and 3cr (99,73%) credible regions have PDF levels of 0.53, 0.10 and 0.01, 
respectively. As Q^o + ^m0 + &c0 = 1, the probability for a spatially flat Universe is on the line Q m o + f^cO = 1> above it 
we have the region for a closed Universe (k > 0, UkO < 0), and below, the region for an open Universe (fc < 0, QkO > 0)- 
The integral of the joint PDF over the upper region gives the likelihood of a closed Universe being 84.0%, althought the 
joint PDF peak is located in the open Universe region. 



way, the case of (Qko, ^mo) = (0, 0) is favoured at a confidence level of 99.3% (2.69cr). 



5 Conclusion 

In the present paper constraints on a Chaplygin gas model using type la supernovae data were settled out 
using a Bayesian statistics. The model contains as matter content the Chaplygin gas and a pressureless 
fluid, whose density parameters are Sl c o and Q m0 respectively. The curvature term Qko, the sound velocity 
of the Chaplygin gas A and the Hubble parameter H are also free parameters. At the 2<t level, the 
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results indicate H = 62. ltf l km /Mpc.s, Q k0 = -OMt\H, « m o = O.Otgjf , fl cQ = 1AQ±{[\1, A = c 2 s = 
0. 93t^2i c i *o = 14-2li;|Gt/ and q = — 0.98l l62- This results must be compared with those obtained 
replacing the Chaplygin gas by the cosmological constant (A = 1), for which Ho = 62.2 ± 3.1 km/Mpc.s, 
n k0 = -O.SOtlil Vmo = 0.58t°;tf , fi A = 1.21±°;|i, t = 15A±H Gy and q = -0.991°;™ 

The CGM is more likely than the ACDM model with a confidence level of 55.3% when all free 
parameters are considered. This percentage is much higher if some restriction on the parameters are 
imposed: fixing fl k Q = 0, the CGM becomes preferred with respect to ACDM at 63.5%; this confidence 
level mounts to 91.8% if, instead, fi m o = 0. The Bayesian analysis of the flat case with Chaplygin gas 
and baryonic matter (f2bo = 0.04) implies that the CGM is preferred with respect to the ACDM model 
at 99.4% confidence level. In general, the predicted value for the sound velocity of the Chaplygin gas is 
close to the cosmological constant value (A = 1), but after integrating on the various parameters the 
ACDM case becomes quite disfavoured. 

Other important differences between the CGM and ACDM concern the pressureless matter parameter 
and the age of the Universe. For the former, the CGM favours a zero value for this cold dark matter 
component, in agreement with the idea that the Chaplygin gas may unify dark matter and dark energy, 
and in contrast with ACDM where a non-negligible fraction of the matter in the Universe must appear 
under the form of dark matter. 

A very remarkable discrimination between the CGM and ACDM occurs when the pressureless matter 
parameter represents the baryonic matter, fl^o = 0.04 and the Universe is spatially flat. This case is 
favoured at 98.6% of confidence level for the CGM, while for the ACDM model it is excluded with 
99.6% of confidence level. This result renders the CGM quite attractive in view of the predictions 
of almost all primordial inflationary scenarios, which lead to a flat Universe, and also in view of the 
unification program for dark energy and dark matter through the Chaplygin gas. 

One of the most important conclusions of this work is that the predicted value for the dark matter 
parameter Q m o is peaked in the zero value. This reinforces the idea that the Chaplygin gas may unify 
dark matter and dark energy as its behaviour in terms of the scale factor suggests. This unification 
program has been criticized |20U21| because, besides some other reasons, the matter power spectrum in a 
pure CGM exhibits oscillations that are not observed in the recent 2dFGRS data, and the phase space 
of possible configuration is highly concentrated around the cosmological constant value. This seems to 
be a strong argument against the unified model. However, we must remark that the authors of Ref. |2()j 
employ an one fluid model. Even if the CGM may unify dark matter and dark energy, baryons exists 
anyway, even if in a small fraction (£lbo ~ 0.04). This may seems to be irrelevant, but the behaviour 
of a two fluid models is generally very different from an one fluid model (see, for example, Ref. |22|'). 
and this point deserves, in our opinion, a deeper analysis. In what concerns the results of Ref. |21| . we 
observe that the crossing of all observational data, including type la supernovae, may put the Chaplygin 
gas again at a competitive level with respect to ACDM . Our results indicate clearly that the Chaplygin 
gas is preferred with respect to ACDM if only type la supernovae data is taken into account. 

The CGM predicts a Universe younger than ACDM , but still in agreement with other astronomical 
data, in particular with the age of globular clusters. Both models favour a closed Universe, with a very 
small difference for the value of the curvature parameter. On the other hand, the value of the Hubble 
parameter is essentially the same in both models, as well as the deceleration parameter qo, for which the 
data indicate a highly negative value, near —1. The value of the deceleration parameter becomes less 
negative for a flat Universe and when CDM is absent. 

In Ref. |38| , statistics of gravitational lenses where used to constrain the proportion of dark matter 
and Chaplygin gas in a flat Universe, and the authors found fi TO o ~ 0.2. We remark also that our results 
are consistent with the particular case described in Ref. |19| . 

These results may be compared with those recently obtained from the WMAP observatory for the 
spectrum of the anisotropy of the cosmic microwave background radiation. Using also the 2dFGRS 
and Lyman-a forest data, and fixing a ACDM model, it has been obtained that f2 m o/i 2 = 0.133 ± 
0.006, h = 0.72 ± 0.03, where H Q = lOOhMpc/km.s, and n k0 = 0.02 ± 0.02. We notice that the 
dark matter component is much smaller than the value deduced from the supernovae data for ACDM 
(!] m o = 0.58lggg), and at same time the Hubble parameter is greater. The smaller uncertainty on the 
dark matter component is natural in this case, since the CMB spectrum gives a better estimation on the 
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total matter in the Universe through the position of the first acoustic peak, and it is verified that an 
almost flat Universe is favoured by the data, for the ACDM model, in contrast with what the supernovae 
data indicate for the same theoretical framework. 

For the CGM there is a large uncertainty also in the estimation of the curvature parameter if only 
supernovae data are used. But, its quite remarkable that, if the dark matter density is fixed to zero in 
this model the type la supernovae data indicates an almost flat Universe: flko = 0.17^5 si- ^ must be 
stressed that a proper comparison between the WMAP results with our analysis for the Chaplygin gas 
requires that the WMAP data must be analysed using the CGM, what is one of the natural extension 
of the present work. 

We hope that, in the future, more SNe la data (from the SNAP project [351001! etc) with small 
observational errors will impose stringent constraints on the parameter estimation for the CGM and 
the ACDM model, so the parameter credible regions become narrow enough to rule out one of these 
cosmological models. For example, the estimation of the parameter A could favour one of the models with 
high confidence level (> 2a), the estimation of £l m o and to could be incompatible with other independent 
and well accepted estimations therefore excluding some cosmological models, etc. 
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